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Detailed  statistical  rate  calculations  combined  with  electron  capture  theory  and  kinetic  modeling  for 
the  electron  attachment  to  SF6  and  detachment  from  SF6"  [Troe  et  al. ,  J  Chem.  Phys.  127,  244303 
(2007)]  are  used  to  test  thermionic  electron  emission  models.  A  new  method  to  calculate  the  specific 
detachment  rate  constants  kdc[(E)  and  the  electron  energy  distributions  /(E,e)  as  functions  of  the 
total  energy  E  of  the  anion  and  the  energy  e  of  the  emitted  electrons  is  presented,  which  is 
computationally  simple  but  neglects  line  structures  in  the  detailed  kdcl(E).  Reduced  electron  energy 
distributions  f(E,E/(e))  were  found  to  be  of  the  form  (e/(e))M  exp(-£/(£))  with  /i~ 0.15,  whose 
shape  corresponds  to  thermal  distributions  only  to  a  limited  extent.  In  contrast,  the  average  energies 
(e(E))  can  be  roughly  estimated  within  thermionic  emission  and  finite  heat  bath  concepts.  An 
effective  temperature  Td(E)  is  determined  from  the  relation  £'-EA  =  (£'Sp6(r^))'f kTd,  where 
<E$P6(rf/))  denotes  the  thermal  internal  energy  of  the  detachment  product  SF6  at  the  temperature  Td 
and  EA  is  the  electron  aflinity  of  SF6.  The  average  electron  energy  is  then  approximately  given  by 
( s(E))=kTd(E ),  but  dynamical  details  of  the  process  are  not  accounted  for  by  this  approach. 
Simplified  representations  of  fcde,(f)  in  terms  of  Td(E )  from  the  literature  are  shown  to  lead  to  only 
semiquantitative  agreement  with  the  equally  simple  but  more  accurate  calculations  presented  here 
An  effective  “isokinetic”  electron  emission  temperature  Te(E)  does  not  appear  to  be  useful  for  the 
electron  detachment  system  considered  because  it  neither  provides  advantages  over  a  representation 
of  kdc[(E)  as  a  function  of  Td(E),  nor  are  recommended  relations  between  Te(E)  and  Td(E)  of 
sufficient  accuracy.  ©  2009  American  Institute  of  Physics.  [DOF  10.1063/1 .3 149782] 


I.  INTRODUCTION 

Molecular  systems  with  several  internal  degrees  of  free¬ 
dom  often  act  as  their  own  heat  bath.  In  this  picture,  the  lixed 
microcanonical  internal  energy  E  is  distributed  over  the  in¬ 
ternal  degrees  of  freedom  in  a  quasithermal  manner  and  of¬ 
ten  it  is  expedient  to  deline  an  effective  temperature  T 
through  the  relation 

< E(T))  =  E ,  (1.1) 

( E{T ))  here  denotes  the  average  thermal  internal  energy  of 
the  molecules  in  a  canonical  ensemble  at  the  temperature  T. 
In  other  words,  an  effective  temperature  T  is  defined  such 
that  the  average  energy  ( E(T ))  in  a  canonical  ensemble 
equals  the  microcanonical  energy  E. 

Equation  (1.1)  allows  one  to  compare  physical  properties 
Y  of  the  molecular  system  in  microcanonical  and  canonical 
ensembles.  Some  properties  depend  on  the  internal  energy 
distribution  in  such  a  way  that  YiE)  and  Y{T)  approach  each 
other  provided  that  a  valid  relation  (1.1)  is  fulfilled.  For  in¬ 
stance.  it  has  been  demonstrated  both  experimentally1  and 
theoretically^ 3  that  UV  absorption  cross  sections  rr(\)  of  vi- 
brationally  excited  polyatomic  molecules  in  microcanonical 
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and  canonical  ensembles  are  practically  the  same  when 
(ri\,E)  and  (r(\,T)  are  compared  and  Eq.  (1.1)  is  obeyed. 
The  situation  may  change  when  product  energy  distributions 
of  bond  fission  processes  are  considered.  As  a  consequence 
of  the  heat  bath  concept,  one  would  expect  quasithermal  dis¬ 
tributions  of  product  translations,  rotations,  and  vibrations 
which  all  correspond  to  the  same  temperature,  see,  e.g  Ref. 
4.  However,  dynamical  effects  often  are  found  to  produce 
differences  between  the  respective  effective  temperatures, 
see,  e.g.,  the  calculations  for  the  fragmentation  of  n-propyl 
benzene  cations  in  Ref.  5.  The  concept  finally  loses  its  sig- 
nilicance  when  the  physical  property  Y  samples  the  high- 
energy  tail  of  the  internal  energy  distribution.  In  canonical 
ensembles  the  energies  extend  to  high  values,  whereas  they 
are  restrained  by  the  value  of  E  in  microcanonical  en¬ 
sembles.  Correspondingly,  specilic  rate  constants  k(E)  and 
their  thermal  analogs  k( 73  of  processes  with  large  threshold 
energies,  in  general,  cannot  be  related  by  means  of  Eq.  (1.1). 

The  present  article  deals  with  a  phenomenon  which  has 
been  also  analyzed  in  terms  of  the  heat  bath  concept,  i.e.,  the 
detachment  of  electrons  from  metastable  molecular  anions  or 
negatively  charged  clusters.  The  phenomenon  has  been 
termed  “thermionic  electron  emission;”  for  an  extensive  re¬ 
view  see,  e.g..  Ref.  6.  Once  again  the  question  arises  as  to 


0021  *9606/2009/1 30<  24 )/244303/12/$25.00  130,  244303-1  ©  2009  American  Institute  of  Physics 

Downloaded  29  Jun  2009  to  146.1 53.144  35.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jcp.aip.org/jcp/copyright.jsp 


244303-2  Troe,  Miller,  and  Viggiano 


J  Chem.  Phys  130,  244303  (2009) 


how  far  the  heat  bath  concept  holds  or  to  what  extent  dy¬ 
namical  effects  require  modifications.  This  is  the  subject  of 
the  present  article. 

We  base  our  analysis  of  electron  detachment  (or  “emis¬ 
sion”)  rates  on  a  theoretical  treatment  of  the  detachment  dy¬ 
namics  by  statistical  ummolecular  rate  theory  combined  with 
electron  capture  theory.  In  addition,  specific  kinetic  effects  of 
the  considered  system,  which  influence  the  attachment/ 
detachment  process,  are  accounted  for.  An  essential  element 
of  the  treatment  is  the  link  between  the  lifetimes  of  the  ex¬ 
cited  anions  with  respect  to  electron  detachment  and  the 
cross  sections  for  the  corresponding  electron  attachment. 
This  link  is  provided  by  detailed  balancing.  It  is  automati¬ 
cally  included  in  statistical  rate  theory  and  it  is  common  to 
most  thermionic  emission  models,  see,  e.g.,  Ref  6.  A  neces¬ 
sary  condition  for  the  application  of  detailed  balancing,  how¬ 
ever,  is  that  the  linked  processes  are  truly  reverse.  This  as¬ 
pect  will  also  be  examined  in  the  present  article. 

Within  a  finite  heat  bath  concept,  effective  temperatures 
of  the  parent  anions,  the  daughter  neutrals  and  the  emitted 
electrons  have  been  formulated  in  a  number  of  ways.6-10 
These,  however,  are  not  necessarily  consistent  with  Eq.  (1  1). 
As  the  characterization  of  electron  attachment  and  detach¬ 
ment  has  made  progress  over  the  last  years,  it  appears  useful 
to  reconsider  these  effective  temperatures  and  to  investigate 
the  influence  of  dynamic  factors. 

It  is  clear  that  the  definition  of  an  effective  electron 
emission  temperature  at  low  electron  energies  meets  with 
problems.  For  illustration  we  consider  a  process  of  the  type 

A-*  — >  A  +  (1.2) 

at  an  energy  E  of  the  excited  parent  anion  A“*  which  is  not 
much  larger  than  the  threshold  energy  of  the  process.  The 
latter  is  assumed  to  be  given  by  the  electron  afhnity  (EA)  of 
A.  If  the  energy  E-EA  available  to  the  electron  is  smaller 
than  the  energy  of  the  first  vibrationally  excited  state  of  the 
daughter  neutral  A,  there  is  no  multilevel  heat  bath  of  A  over 
which  the  energy  E-EA  could  be  distributed.  Instead  the 
electron  energy  e  is  fixed  and  given  by  the  difference  f,  =  E 
-EA,  whereas  the  internal  energy  of  A,  E-EA-e,  is  zero. 
Although  Eq.  (1.1)  could  be  applied  to  this  situation  as  well, 
the  derived  temperature  would  have  only  little  in  common 
with  that  of  a  heat  bath.  On  the  other  hand,  if  E-EA  reaches 
up  into  the  vibrational  quasicontinuum  of  A,  there  will  be 
broad  distributions  of  electron  energies,  finite  heat  bath  mod¬ 
els  might  be  considered,  and  an  effective  electron  tempera¬ 
ture  could  be  defined  more  meaningfully.  Two  questions 
arise;  At  which  energies  E-EA  does  the  transition  from  low- 
to  high-energy  behavior  occur,  and  how  are  effective  tem¬ 
peratures  related  to  the  kinetic  and  dynamic  details  of  the 
process  (1.2).  We  focus  attention  on  these  questions  in  the 
present  article  for  a  number  of  reasons.  First,  we  feel  that  the 
statistical  theory  of  electron  attachment/detachment  pro¬ 
cesses  in  the  framework  of  thermionic  emission  models6  of¬ 
ten  has  been  used  in  an  oversimplified  manner.  Second,  we 
recently  have  extended*  Vogt-Wannier  (VW)  and 
Fabnkant-Hotop-type*4  capture  theory  for  electron  attach¬ 
ment.  The  results  of  this  work  may  be  implemented  into 
statistical  calculations  of  electron  detachment  rates.  Third, 


we  have  provided  a  detailed  kinetic  modelling  ,fl  of  the 
mechanism  of  low-energy  electron  attachment  to  SF6  and  the 
corresponding  electron  detachment  from  SF6”, 

SF6“^SF6  +  E\  (1.3) 

This  model  is  “user  friendly,”  it  can  easily  be  generalized  to 
other  attachment/detachment  systems,  and  it  can  incorporate 
specific  dynamical  and  kinetic  factors  for  individual  reac¬ 
tions.  For  reaction  (1.3),  e.g.,  it  was  suggested  that  the  ex¬ 
perimental  cross  sections  for  nondissociative  electron  attach¬ 
ment  include  kinetic  factors  which  account  for  electron 
capture  into  a  virtual  state  e~SFb,  “intramolecular  vibrational 
redistribution”  (IVR)  (or  electron-phonon  coupling)  leading 
from  ^“SF6  to  anionic  SF6  *,  and  inelastic  vibrational  exci¬ 
tation  (VEX)  of  SF6  in  competition  with  IVR  from  <TSF6  to 
SF6  .  VEX  can  either  proceed  through  highly  excited  an¬ 
ionic  SF6  *  or  it  can  avoid  the  intermediate  formation  of 
metastable  SF6  *.  Alternative,  yet  more  complicated,  de¬ 
scriptions  of  the  attachment  process  in  terms  of  Gauyacq- 
Herzenberg  resonance  theory  4  have  also  been  elabo¬ 
rated  and  might  be  compared  with  VW-type  electron  capture 
theory.1114  Finally,  our  recent  measurements  of  thermal  at¬ 
tachment  and  detachment  rates1  through  a  third  law  analysis 
of  the  equilibrium  constant  led  to  a  revised  value  of  the  EA 
of  SF6,  suggested  to  be  EA(SF6)  =  1 ,20(  ±0.05)  eV.  This,  in 
turn,  gave  rise  to  a  revision  of  statistical  calculations  of  the 
specific  rate  constants  kdcl(E)  for  electron  detachment  from 
SF6  *  and  of  the  lifetimes  \/kdcl(E)  of  SF6  *. 

Because  of  the  wealth  of  theoretical  and  experimental 
studies  available  for  the  SF6"  system,  it  appears  attractive  to 
inspect  thermionic  electron  emission  models  for  reaction 
(1.3).  This  is  even  more  timely  since  storage  ring 
measurements”  of  SF.  lifetimes  have  recently  been  evalu- 

°  n  1 

ated  in  terms  of  such  models,”  the  latter  analysis  suggesting 
a  value  of  EA(SF6)=1.4  eV,  which  is  higher  than  the  value 
of  EA=1.20eV  derived  in  Ref.  19.  The  present  article  in¬ 
spects  electron  detachment  rates  from  SF6“  over  a  wide  en¬ 
ergy  range,  employing  electron  capture  theory  combined 
with  statistical  rate  theory  without  the  simplifications  usually 
made  in  thermionic  emission  models.6  It  also  traces  the  rea¬ 
sons  for  the  discrepancy  between  the  EA(SF6)  values  sug¬ 
gested  in  Refs.  19  and  21.  Attention  in  our  work  is  finally 
paid  to  the  transition  from  low-  to  high-energy  behavior. 
With  the  experience  gained  for  the  SF6"  system,  the  treat¬ 
ment  may  be  extended  to  larger  species  such  as  C60“  or  other 
negatively  charged  atomic  clusters.  This  continuation  of  our 
work  will  be  presented  in  future  publications. 

Before  elaborating  the  theory  of  electron  energy  distri¬ 
butions  and  detachment  rates,  one  particular  remark  about 
the  relation  between  cross  sections  for  nondissociative  elec¬ 
tron  attachment  to  SF6  and  detachment  rate  constants  should 
be  made.  It  was  suggested  in  Ref.  15  that  the  experimental 
attachment  cross  sections  to  SF6  contain  contributions  from 
several  kinetic  processes  such  as  primary  electron  capture, 
IVR,  and  VEX.  The  question  arises  which  of  these  contribu¬ 
tions,  through  detailed  balancing,  should  be  included  in  cal¬ 
culations  of  detachment  rate  constants  kdel{E).  Therefore, 
kdct(E)  in  Ref  15  was  calculated  alternatively  with  and  with- 
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out  the  individual  kinetic  factors.  This  is  an  important  kinetic 
aspect  of  the  problem  which  should  be  taken  into  account  in 
all  systems  for  which  thermionic  emission  models  are  ap¬ 
plied.  We  will  readdress  this  point  later. 

II.  STATISTICAL  RATE  THEORY 

In  the  following  we  briefly  summarize  the  basic  equa¬ 
tions  of  statistical  unimolecular  rate  theory1  such  as  they 
were  adapted  in  our  recent  work  to  electron  detachment/ 
attachment  processes.  3  The  starting  point  is  the  specific  rate 
constant  for  detachment, 

k,c[(EJ)=Wde{(EJ)fhpp(EJ)  (2.1) 

at  the  energy  E  and  total  angular  momentum  7,  with  Planck’s 
constant  //,  the  “number  of  open  reaction  channels”  or  the 
“cumulative  reaction  probability”  Wdel(E,7),  and  the  rovibra- 
tional  density  of  states  of  the  parent  anion  pp(EJ).  One  may, 
or  may  not,  include  in  VVdel(£,7)  the  electronic  degeneracy 
factor  of  2,  originating  from  the  spin  of  the  emitted  electron 
because  the  same  factor,  now'  originating  from  the  doublet  of 
the  anion,  is  contained  in  pp(EJ).  These  factors  cancel  in 
Eq.  (2.1).  It  appears  confusing  to  explicitly  account  for  this 
factor  in  lVdel(£,7)  and  to  leave  it  implicitly  in  pp(EJ). 
Since  they  cancel,  in  the  present  work  the  electronic  degen¬ 
eracies  are  omitted  from  the  beginning  both  from  Wdel(£,7) 
and  (>p(E,J) 

Wdel(£,7)  is  given  by 

wde, ( E,J)  =  £/>(£  (2.2) 

where  the  £(),  are  the  energy  levels  of  the  daughter  neutral 
species  and  E(£-£0,)  are  the  transmission  coefficients  for 
the  outgoing  (in  detachment)  or  incoming  (in  attachment) 
electrons.  Because  of  microscopic  reversibility,  the  latter  are 
the  same  under  the  condition  that  the  same  dynamical  phe¬ 
nomena  are  compared  The  cross  sections  for  attachment, 
starting  from  the  level  £()/<  then  are  directly  related  to  the 
detachment  rate  constants  through  detailed  balancing,  i.e., 
through  the  relationship 

X 

<rt*,-)  =  (*/*J,)2(2/+  I  )Pi(kpl).  (2.3) 

/=0 

where  kpi-pj fi  =  \2jx(E-E{)l)l ti  and 

/»(£-£«)  =  2(2/  +  I  )Pl(kpt).  (2.4) 

/=0 

Above,  kpi  is  the  wavenumber  for  momentum  p,  and  fi  is  the 
reduced  mass  of  SF6  and  e~.  If  only  j-waves  need  to  be 
considered,  such  as  in  the  electron  detachment/attachment  of 
SF6",  one  has  P{E-E0i)  =  Pl=i)(kpi).  This  property  of  the  SF6" 
system  very  much  simplifies  the  treatment.  If  higher  waves 
are  included,  their  angular  momentum  (quantum  number  /) 
contributes  to  the  total  angular  momentum  (quantum  number 
7)  which  needs  particular  attention  when  Eqs.  (2.3)  and  (2.4) 
are  combined  with  Eq.  (2.2). 

Electron  capture  or  resonance  theories  1  s  provide  suit¬ 
able  expressions  for  the  P{(kpt),  both  for  the  case  of  daughter 


neutrals  having  induced  or  having  permanent  dipole 
moments.1  The  study  of  the  SF6~  system  in  Ref.  15  sug¬ 
gested  that  the  mechanistic  complications  arising  from  IVR 
and  inelastic  VEX  can  empirically  be  accounted  for  by 
factors  to  be  multiplied  with  the  Pjik^).  Below  we  inspect 
the  consequences  of  such  effects  for  electron  emission 
properties. 

For  each  reaction  channel  /,  the  energy  of  the  electrons 
E;  is  equal  to  the  energy  difference  E~E0h  i.e., 

e,  =  E-  E0l.  (2.5) 

In  the  following,  we  are  interested  in  partial  specific  rate 
constants  £deJ(£,c,)  for  electron  detachment.  As  long  as  only 
one  channel  E0i  is  open,  i.e.,  £0,  =  EA,  one  has  e,=8  =  E 
—  EA  and  kdel(E,F,t)  is  given  by 

*de[(£,r,)  =  *dcl(£»  =  Pie)lhpp(E).  (2.6) 

[In  the  following  we  omit  the  7-dependence  which  was 
elaborated  in  Ref.  15  and  only  consider  A:det(E,7=0)  which 
we  denote  by  &dcl(£).  However,  we  note  that,  for  .?-wave 
scattering,  one  has  7=/,  where  j  corresponds  to  the  total 
angular  momentum  of  the  neutral  daughter.  In  principle,  all 
kdc:(EJ)  need  to  be  considered,  see  Ref.  15].  In  the  low- 
energy  region  there  is  no  difference  between  partial, 
/cdel(£,s),  and  total,  £de,(E),  detachment  rate  constants.  In¬ 
creasing  the  energy  into  the  sparse  manifold  of  the  daughter 
vibrational  levels  £0(,  at  a  given  total  energy  £  there  are 
several  levels  /  with  electron  energies  e{  which  contribute, 
and  one  has 

*de,(£)  =  2*de.<i».  (2-7) 


with 


*det(£,e,)  =  £(e,)//tp,,(£) . 


(2.8) 


In  this  case  one  will  have  a  discrete  distribution  F(£,e()  of 
electron  energies  defined  by 


F(E,e,)  =  W£.b,.)/W£-)  =  £(«,■) 


2>(e,). 


(2.9) 


Further  increasing  the  energy  into  the  vibrational  quasicon¬ 
tinuum  of  the  daughter  states,  one  switches  to  a  quasicon- 
tinuous  distribution  function  /(E,e)  of  electron  detachment 
rate  constants  defined  by 

/t£,8)sW£,fi)/U£).  (2.10) 

In  this  case,  the  summation  over  i  in  kdc[(E)  is  replaced  by  an 
integral  over  e  weighted  by  the  vibrational  density  of  states 
Pd{E- EA-e)  of  the  daughter  neutral.  The  combination  of 
Eqs.  (2. 1 )— ( 2.5)  then  leads  to  the  partial  rate  constant  per 
energy  e  interval. 


kdJE,e) 


fipd( E  -  EA  -  e)cr(£)E 
7 rt>'pp(E) 


(2.11) 


and,  hence,  to  the  distribution  function  per  energy  interval. 
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aU]epjE-EA-e) 

1  ,£  ~  fE0~E\r(r,)*pd(E-EA-r.)dF/ 


(2.12) 


Equation  (2.11)  agrees  with  the  general  formulation  in  Ref. 
6  However,  one  should  note  again  that  in  the  present  work 
pd  and  pp  both  are  pure  vibrational  densities  of  states  and  the 
cancelling  electronic  degeneracies  have  been  taken  care  of 
separately.  [The  factor  of  2  in  the  numerator  of  Eq.  (2)  of 
Ref.  6  then  disappears,  see  above.]  In  addition, 

/(£,e),  and  er(e)  should  correspond  to  the  same  dynamical 
and  kinetic  process  in  the  forward  and  reverse  directions. 
This  caveat  is  important  to  be  kept  in  mind 

In  the  following,  we  explicitly  determine  *<.  «(£), 
£(£,£*),  and  /(£,£)  for  electron  detachment  from  SF6"  and, 
we  analyze  the  contributions  of  the  various  kinetic  factors  to 
these  quantities.  We  base  our  calculations  on  the  molecular 
parameters  of  SF6"  and  SF6  from  Ref.  22,  such  as  given  in 
Appendix  and  also  employed  in  Ref.  15. 


III.  DETACHMENT  RATES  AND  ELECTRON  ENERGY 
DISTRIBUTIONS  FROM  SF6" 

Total  specific  rate  constants  A:det(£,7)  for  electron  de¬ 
tachment  from  SF6“  in  Ref.  15  were  calculated  in  detail  and 
documented  over  the  range  £=1.2-2.2eV.  The  results 
showed  a  considerable  fine  structure  caused  by  the  sparse 
character  of  the  vibrational  manifold  of  the  £0/  of  SF6.  The 
7-dependence  was  also  illustrated  over  the  range  for  J 
=0-200.  Calculations  were  compared  employing  (i)  P(E 
-£0,)  =  1  at  E 5* E0l  [and  P(E-E0i)  =  0  at  £<£0J,  00  £(£ 
-Eoi)  from  pure  electron  capture  theory,  and  (iii)  P(E-E0i) 
including  kinetic  factors  arising  from  IVR  and  VEX.  In  the 
following,  we  analyze  the  corresponding  kdel(E)  and  the  en¬ 
ergy  distribution  of  the  emitted  electrons,  alternatively  in¬ 
cluding  or  omitting  the  kinetic  factors. 

In  our  analysis  we  follow  the  approach  outlined  in  Ref. 
15.  We  start  with  the  “Rice-Ramsperger-Kassel-Marcus 
(RRKM)-type”  treatment  where  the  P(E-E0l)  are  step  func¬ 
tions  rising  from  zero  to  unity  at  all  individual  £2*£0i.  We 
then  continue  with  pure  electron  capture  theory,  improving 
Klots’  analytical  approximation  to  P(E-  E0l)  for  s-wave 
capture  within  the  VW  model,  using  the  expression  derived 
in  Ref.  1 2, 

P{E- Eoi)  =  Pvw(£  -  E0i) 

555  1  -  0.5  exp(-  2ki)  -  0.5  exp(-  6/q),  (3.1) 

with 

Kj  =  fie\2 a{E  -  E0i)/h~ .  (3.2) 

Here,  e  notes  the  electronic  charge  and  a  is  the  isotropic 
polarizability  of  SF6,  taken  as  a=6.54  X  10"24  cm3.  This 
leads  to  the  reduced  energy  scale  #q=  1 .80[(£-£0f)/eV]1  2 
Suggesting1  that  both  the  experimental  nondissociative  at¬ 
tachment  cross  sections"4-2^  and  the  thermal  detachment  rate 
constants15,2  2  for  SF6~  are  not  only  determined  by  pure 
electron  capture  but  also  include  contributions  from  IVR,  we 
then  modify  £(**,)  by  an  additional,  empirical,  factor 
P1vr(k,).  In  Ref.  15  this  factor  was  assumed  to  have  the  form 


PIVR(/c,)  =  exp(-c,*f).  (3.3) 

The  analysis  of  experimental  thermal  attachment  rate  con¬ 
stants  between  50  and  700  K  led  to  a  temperature  dependent 
fit  parameter,  being  1.92  for  300  K.  While  the  parameter 
C/  is  still  quite  uncertain  for  the  range  of  50-200  K,  its  value 
appears  well  established  for  200-700  K,  see  the  analysis  in 
Ref.  15.  It  should  nevertheless  be  emphasized  that  the  kinetic 
factor  of  Eq.  (3.3),  as  suggested  by  the  analysis  of  thermal 
attachment  rate  constants,  is  based  on  information  for  a  lim¬ 
ited  energy  range  only,  One,  therefore,  cannot  expect  that  Eq. 
(3.3)  applies  to  larger  energy  ranges.  We  keep  this  restriction 
in  mind  w  hen  larger  energy  ranges  are  considered.  [Equation 
(3.3)  may  also  include  contributions  from  an  energy  depen¬ 
dence  of  the  effective  polarizability  of  SF6  which  is  not  fur¬ 
ther  considered  here.] 

The  analysis  of  experimental  nondissociative  attachment 
cross  sections  suggested  that  yet  another  dynamical  factor 
should  be  included  in  P(Kj)  for  attachment  which  accounts 
for  inelastic  collisions  between  e~  and  SF6  without  SF6~* 
being  formed.  This  factor  was  suggested  5  to  be  of  the  form 
PVEX=1  for  k,^k^  and 

/>VEX(*/)«expl-c2(K?-K^)J  (3.4) 

for  where  =0.588  corresponds  to  an  energy  E 

-EA  of  0  096  cV  and  c2  empirically  was  fitted  as  c2^6. 0 
for  300  K.  If  this  assumption  of  competing  electron  attach¬ 
ment  forming  SF6  *  and  inelastic  electron  scattering  avoid¬ 
ing  SF6  *  is  correct,  then  pVEX(^.)  should  be  included  in  the 
overall  nondissociative  attachment  cross  sections,  but  it 
should  not  be  included  in  the  detachment  rate  constants.  On 
the  other  hand,  P^r(k,)  should  be  included  in  the  detach¬ 
ment  rate  constants  because  detachment  also  invokes  IVR. 
As  this  interpretation,  however,  is  still  uncertain,  /cdet(£)  in 
Ref.  15  alternatively  was  calculated  with  and  without  the 
factor  PVEX(Aq)  In  the  following,  we  omit  the  PVEX(/c<)  fac¬ 
tor,  assuming  that  the  described  mechanism  is  correct.  We 
emphasize  that  experimental  nondissociative  attachment 
cross  sections  and  detachment  rate  constants  then  are  not 
directly  linked  by  detailed  balancing  because  they  only 
partly  correspond  to  reverse  kinetics.  [In  Refs.  15  and  16 
/>,VRU,)  and  PVEXU,)  could  be  separated  only  up  to  an  en¬ 
ergy  of  about  E-Eq^0.2  eV,  while  for  larger  energies  up  to 
about  £-£0~  2  eV,  only  the  product  of  the  two  factors  was 
fitted  to  the  experimental  total  attachment  cross  sections;  a 
separation  of  the  two  factors  at  large  energies  was  not  pos¬ 
sible.  This  limits  the  possibilities  to  design  /cdet(£).  It  is  even 
not  proven  that  the  differences  between  experimentally  ob¬ 
served  thermal  nondissociative  attachment  rate  coefficients 
at  temperatures  below  about  700  K  and  calculations  with  the 
VW  electron  capture  model  are  due  to  IVR  or  as  well  to 
VEX.  In  the  latter  case,  kdcl(E)  would  only  contain  VW  con¬ 
tributions.]  In  spite  of  the  remaining  interpretational  uncer¬ 
tainties,  Eqs.  (2.3),  (2.4),  and  (3. 1 )— (3.4)  in  Refs.  15  and  16 
were  shown  to  provide  a  very  good  representation  of  the 
available  experimental  attachment  cross  sections  and  thermal 
attachment  rate  coefficients  [both  ka[(T^as=Te])  and  kal(Tg2iS 
*7e|)]. 
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FIG.  I  Discrete  distributions  F(£,e)  of  electron  energies  e  in  the  detachment  from  SF6~  at  the  energies  £— EA=0.025  eV  (a).  0.05  eV  (b).  and  0.1  eV  [(c) 
and  (d )].  (d)  shows  A F(E.tt)  Ak  with  bins  A*:  =  0.02  eV:  the  continuous  curve  is  from  Fq.  (2.12).  see  text 


With  the  given  expressions,  specific  rate  constants 
kdel(E)  as  well  as  the  distributions  £(£,£,)  and  /(£,e)  of  the 
electron  energies  defined  in  EUqs.  (2.9)  and  (2.10)  can  directly 
be  calculated.  First,  at  energies  E-EA=s  0.0429  eV  (cor¬ 
responding  to  the  first  excited  vibrational  state  of  SF6  at 
£()f  =  0.0429  eV=/rc  346  cm-1),  one  has 


F(e)  =  1 


(3.5) 


pAE') 


[T  +  a(£)£j2! 

(s-WUhP,  ' 


(3.7) 


where  hvi  are  the  5  vibrational  quanta  and  Ez=  I  /2X-=j hVj  is 
the  vibrational  zero  point  energy  of  $F6.  The  Whitten- 
Rabinovitch  correction  function  a(E')  is  given  by4 

a(E')=  I  (3.8a) 


with 


Equation  (3.5)  corresponds  to  the  range  where  only  a  single 
channel  is  open.  Second,  in  the  range  of  the  sparse  manifold 
of  vibrational  levels  of  SF6,  the  contributions  from  individual 
channels1  E{)l  are  easily  summed  up  individually.  Figure  1 
shows  examples  for  three  values  of  E-EA  in  the  range  of 
0-0. 1  eV.  As  PVEX(  Kj)=  1  over  the  range  of  0-0.096  eV,  the 
question  of  including  or  excluding  EVEX(/q)  practically  does 
not  arise  here.  Figure  1  well  illustrates  the  discrete  nature  of 
E(£,£j)  which  has  its  equivalent  in  the  marked  fine  structure 
of  the  calculated  Adel(£),  see  below'. 

Moving  up  to  higher  energies,  into  the  vibrational  qua¬ 
sicontinuum  of  SF6,  the  distribution  function  /(£,e)  of  Eq. 
(2.10)  with  Eq  (2.12)  takes  the  form 


fiE.s) 


P(E,e)pj(E-  EA-e) 
/TF'V(£.e)P,,(£  -  EA  -  e)de ' 


(3.6) 


The  use  of  Eq.  (3.6)  requires  a  smoothed  expression  for  the 
density  of  states  p,/(£').  Most  conveniently  this  is  obtained 
from  the  Whitten-Rabinovitch  approximation,4 


log,0  o> «  -  1 .0506(E7£ya2!S  at  E’  >  £.,  (3.8b) 


w 1  «  5 (E'/E.)  +  2.73(E72.)0-5  +  3.5 1  at  E'  <  E, 

(3.8c) 


and 


P  =  (s-\)2'2(hv,)2/s 

/=  I 


(3.8d) 


At  the  high-energy  end  of  e,  i.e.,  at  e— >E-EA,  p,/(£') 
reaches  down  into  the  sparse  manifold  of  vibrational  energy 
levels  of  SF6  However,  f(e)  markedly  decreases  with  e  at 
the  high-energy  end  such  that  this  is  not  a  practical  problem 
A  smooth  transition  between  the  continuous  f(s)  and  the 
discrete  AF(k,)/ As  is  obtained  when  Eqs.  (3.6),  (3.7),  (3.8a), 
(3.8b),  (3.8c),  and  (3.8d)  are  used.  For  example.  Fig.  1(d) 
sums  the  discrete  contributions  to  F(E,e)  of  Fig.  1(c)  into 
bins  of  width  A£=0.02  eV  and  compares  the  corresponding 
AE(£,£,)/Ae  with  /(E,e),  the  latter  being  calculated  with 
Eqs.  (3.6),  (3.7),  (3.8a),  (3.8b),  (3.8c),  and  (3.8d).  The  ap- 
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FIG.  2.  Specific  rate  constants  k^iE)  for  electron  detachment  from  SF6“.  Curves  from  top  to  bottom:  (i)  RRKM-type  model  without  electron  capture 
contributions,  i.e.,  P(E-E0i)  =  1  for  £>£0,;  (ii)  pure  electron  capture  model,  i.e.,  £(£-£0,)  =  Pvw(£-£0();  (iii)  electron  capture  with  IVR  contribution,  i.e  , 
P(E-£0,)  =  Pvw(£“£0()PIVR(£-£0i);  (iv)  electron  capture  with  IVR  and  VEX  contributions,  i.e.,  P(£-£O()  =  £vw(£-£o/)/>fVR(E-Eo()PVEX(£'-^0/)-  Full 
curves:  Detailed  calculations  from  Ret.  15;  dashed  curves  for  the  upper  three  models  Simplified  calculations  with  F>q.  (3.9)  from  the  present  work,  see  text. 


proach  of  the  discrete  results  by  the  continuum  distribution  is 
well  documented. 

kdct(E)  in  Ref.  15  with  Eqs.  (2.1W2.4)  and  (3.1)-(3.4) 
was  calculated  without  further  simplifications.  With  increas¬ 
ing  number  of  contributing  channels  i  at  increasing  energy, 
however,  the  calculation  becomes  increasingly  cumbersome. 
The  explicit  calculation  then  is  considerably  simplified  when 
the  Whitten-Rabinovitch  approximation  for  the  densities  of 
states  is  employed  and  A:det(^)  is  derived  from  the  expression 

rE-  EA 

P(E,8)pd(E-EA-e)delhpp(E).  (3.9) 
Jo 

However,  one  sacrilices  the  line  structures  of  kdc[{E).  Figure 
2  compares  the  simplified  £det(£)  from  Eq.  (3.9)  with  de¬ 
tailed  results  from  Ref.  15  such  as  calculated  with  discrete 
channels  /.  Apart  from  the  fine  structure  the  agreement  looks 
very  satisfactory.  Employing  Eqs.  (3.7),  (3.8a),  (3.8b),  (3.8c), 
and  (3.8d),  there  is  a  minor  difference  (about  30%)  between 
the  simplified  and  accurate  results  which  reflects  the  inaccu¬ 
racy  of  Eq.  (3.8a),  (3.8b),  (3.8c),  and  (3.8d)  at  small  ener¬ 
gies.  In  part  this  is  due  to  the  degeneracy  of  the  vibrational 
frequencies  used  in  the  calculations  (see  Appendix).  It  is 
known  that  the  Whitten-Rabinovitch  approximation  then  is 
less  accurate  than  that  is  normally  the  case.  One  can  mostly 
remove  this  problem  by  slightly  modifying  the  Whitten- 
Rabinovitch  parameters  p  in  Eq.  (3.8d).  We  have  done  this 
by  fitting  the  Whitten-Rabinovitch  densities  of  states  for  SF6 
near  £-EA  =  0.29(  ±  0.05)  eV  to  Beyer-Swinehart  counting 
results.4  This  is  achieved  by  reducing  the  parameter  p  from 
Eq.  (3.8d)  by  a  factor  of  0.907  65.  [As  E  is  much  larger  than 
£-EA,  pp(E)  needs  no  such  correction.]  The  inaccuracy  of 


Eq.  (3.9)  then  nearly  disappears.  Figure  2  compares  the  cor¬ 
responding  results  with  three  of  the  detailed  calculations 
from  Ref.  15.  Figure  2  shows,  from  top  to  bottom,  &del(£)  for 

(i)  the  “RRKM-type  model”  [i.e.,  P(E-E0i)=  1  for  E^E0il 

(ii)  for  the  pure  electron  capture  “VW-type  model”  [i.e., 
P(E-E0i)=P™(E-E0i )  from  Eq.  (3.1)],  (iii)  for  a  “VW 
+  IVR-type  model”  with  electron  capture  followed  by  IVR 
[i.e.,  P(E-En/)  =  Evw(£-E0,)EtvR(E-E0<)  from  Eqs. 
(3. 1 )— (3.3)],  and  (iv)  for  a  “VW  +  IVR  +  VEX-type  model” 
with  electron  capture  followed  by  IVR  and  competing  with 
VEX  [i.e.,  E(E-E0<)  =  Pvw(E-E0/)  =  P1vr(E-E0<)/>vex(E 
-£0,)  from  Eqs.  (3. 1 )— (3.4)]. 

Having  gained  conlidence  in  the  use  of  the  (slightly  cor¬ 
rected)  Whitten-Rabinovitch  approximation  for  densities  of 
states  in  calculations  of  detachment  rate  constants 
through  Eq.  (3.9),  we  use  this  approximation  also  in  the  cal¬ 
culations  of  electron  energy  distributions  /(£,e)  through  Eq. 
(3.6).  We  first  analyze  the  influence  of  the  various  contribu¬ 
tions  to  E(£,e)  on  the  shape  of  the  distribution  functions. 
We  consider  an  energy  £-EA= 1  cV  which  is  well  above  the 
discrete  range  characterized  in  Fig.  1.  Figure  3  shows  the 
results.  A  representation  of  /(£,e)  as  a  function  of  e  here  is 
chosen.  In  comparing  the  results  with  statistical  unimolecular 
rate  theory4  one  should  remember  that,  unlike  normal  unimo¬ 
lecular  reactions,  the  present  system  is  characterized  by  / 
=  0  such  that  J=  j,  see  above  In  spite  of  this  difference, 
/(£,p)  for  the  RRKM-type  model  looks  similar  to  product 
translational  energy  distributions  for  fragmentations  where 
one  fragment  is  an  atom,  see.  e.g..  Refs.  4  and  29.  The  maxi¬ 
mum  of /(£,e)  then  ts  found  near  e=0.  The  specific  dynami¬ 
cal  effects  of  electron  detachment,  such  as  given  by  electron 
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FIG.  3  Energy  disiribuiion  /(£,*:)  of  electrons  from  ihe  detachment  from  SF6"  (£-EA=l  eV),  full  line  RRKM-type  model,  dashed  line  electron  capture 
VW  model;  dotted  line:  electron  capture  with  IVR  model;  dash-dotted  line  with  peak  near  0. 1  eV:  electron  capture  with  IVR  and  VEX  model,  see  text 


capture  theory,  IVR,  and  VEX  contributions,  change  this  pic¬ 
ture  to  some  extent  such  as  illustrated  in  Fig.  3.  If  the  ex¬ 
perimental  attachment  cross  sections  with  their  marked  drop 
at  £>0.096  eV  would  be  used  for  the  calculation  of  /(E,e), 
then  one  would  have  a  second  maximum  near  e  =  ().  1  eV 
such  as  shown  in  the  figure  However,  as  explained  above  we 
consider  the  curve  corresponding  to  electron  capture  with 
combined  IVR  and  VEX  contributions  as  being  unrealistic, 
because  VEX  presumably  is  not  concerned  in  electron  de¬ 
tachment  Therefore,  in  the  following  we  do  not  further  in¬ 
clude  the  VEX  contributions  from  Eq.  (3.4).  We  also  keep  in 
mind  that  the  IVR  factor  of  Eq.  (3.3)  is  based  on  experimen¬ 
tal  data  over  a  range  of  £  and  E-EA  up  to  about  0.2  eV  only, 
such  that  differences  between  the  VW  and  VW  +  IVR  curves 
at  £>0.2  eV  may  at  least  in  part  be  artificial. 

The  distribution  functions /(E,fi)  of  Fig.  3  can  be  further 
rationalized  by  separately  inspecting  the  average  energy  (fi) 
of  the  emitted  electrons  and  the  shape  of  the  reduced  distri¬ 
bution  functions  /(E, e/(fi)).  We  lirst  illustrate  /(E,fi/(fi)). 
Figure  4  shows  the  curves  of  Fig.  3  in  this  representation. 
Now  the  VW  and  VW+IVR  curves  nearly  coincide  and, 
apart  from  the  maximum  near  £  =  (),  the  RRKM  curve  is  also 
not  too  different  from  the  VW  and  VW  +  IVR  results.  The 
shape  of  the  curves  can  well  be  approximated  by  expressions 
of  the  type 

fi  E,fi/(fi))  *  (£/<£»"  exp(-  fi/(fi')).  (3.10) 

Figure  4  includes  such  curves  for  /i=0.l5,  0.5,  and  I.  Obvi¬ 
ously  the  distribution  with  n  =  0. 15  is  much  closer  to  the 
modelling  results,  both  for  the  VW  and  VW  +  IVR  calcula¬ 
tions,  than  the  “thermal”  distributions  with  n  =  0.5  and  1.  This 
is  confirmed  by  experimental  observations.  For  example,  dis¬ 
tributions  with  n- 0.15,  0.28,  and  0.58  were  observed  for 
electron  emission  from  Cl0”,  C15”,  and  C18"  cluster  anions, 
respectively/  One  may  deduce  from  this  observation  that 
the  distribution  approaches  u  more  thermal  shape  when  the 
cluster  size  increases.  However,  this  effect  may  also  be  due 
to  other  factors  such  as  different  contributions  from  higher 
partial  waves  (/>0).  This  aspect  is  elaborated  in  Part  II  of 
this  work  4  which  applies  the  present  approach  to  larger  car¬ 


bon  clusters  with  varying  shapes.  It  should  be  noted  that  the 
interpretation  of  the  experimental  distributions  in  Ref.  33 
was  done  in  terms  of  orbiting  transition  state/phase  space 
theory,  see  Ref.  4.  This  approach  cannot  be  applied  to  SF6” 
for  at  least  two  reasons:  (i)  for  SF6”,  with  only  s-wave  at¬ 
tachment,  there  are  considerable  angular  momentum  con¬ 
straints,  and  (ii)  dynamical  effects  included  through  electron 
capture  theory  are  neglected.  Nevertheless,  Eq.  (3.10)  pro¬ 
vides  fairly  good  fits  to  the  distribution  functions  /(£,(£)) 
calculated  here,  although  the  differences  between  treatments 
omitting  and  including  electron  capture  theory  (curves  a  and 
b  in  Fig.  4)  are  not  negligible  and  Eq.  (3.10)  does  not  exactly 
reproduce  the  calculated  full  distribution  /(fi,(fi)). 

As  the  various  dynamical  factors  do  not  strongly  influ¬ 
ence  the  shape  of  /(fi,(e)),  it  appears  more  meaningful  to 
focus  attention  on  the  average  energy  (fi).  The  simplicity  of 
Eqs.  (3. 1 )— (3.3),  (3.6),  (3.7),  (3.8a),  (3.8b),  (3.8c),  and  (3.8d) 
allows  one  to  do  thi>  analysis  with  only  little  computational 
effort.  Figure  5  shows  the  dependence  of  (fi)  on  E  and  on  the 
applied  model.  Although  the  general  trend  of  the  three 
curves,  i.e.,  an  increase  in  (fi)  with  increasing  E  is  the  same, 
there  are  considerable  differences  of  (fi)  for  a  given  E.  The 
VW  results  are  always  above  the  RRKM-type  results  The 
VW  +  IVR  results  increasingly  fall  below  the  VW  results. 
However,  as  the  IVR  factor  of  Eq.  (3.3)  was  only  derived 
from  experiments  over  a  range  up  to  about  0.2  eV,  this  may 
at  least,  in  part,  be  due  to  the  limited  energy  range  over 
which  the  IVR  factor  was  derived.  Therefore,  the  VW  curve 
may  be  more  realistic  over  large  energy  ranges.  One  should 
note  again  that  experimental  attachment  cross  sections  at 
higher  energies  differ  from  the  VW  expression  due  to  the 
VEX  contributions  However,  the  latter  should  not  be  in¬ 
cluded  in  detachment  treatments  in  as  far  as  metastable  anion 
states  are  concerned  For  the  curves  a  and  b,  (fi)  is  roughly 
about  one-tenth  of  E.  In  the  following  this  result  is  compared 
with  finite  heat  bath  and  effective  temperature  concepts,  re¬ 
membering  that  (fi)//;  generally  is  assumed  to  be  close  to 
the  effective  daughter  neutral  temperature  T(]  (k  denotes 
Boltzmann’s  constant). 
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FIG.  4  Energy  distnbution  /(£\e(e))  of  electrons  from  the  detachment  from  SF6~  (E-FA=  1  eV).  full  curves=detailed  calculations:  (a)  RRKM-lype  model 
with  (e)=0.0937  eV,  (b)  electron  capture  VW-model  with  (e)=0  109  eV,  and  (c)  VW  +  IVR-model  wtth  (e)=0.0730  eV:  dashed  curves  =  model  functions 
(3.10)  with  exponents  /r  =  0. 1 5  (d),  0.5  (e),  and  1.0  (f). 


IV.  EFFECTIVE  TEMPERATURES 

Our  detailed  calculations  of  detachment  rate  constants 
fcdet(£),  electron  energy  distnbutions  /(£,e),  and  their  aver¬ 
age  ( e(E ))  allow  us  to  analyze  the  appropriateness  of  finite 
heat  bath  concepts  and  the  accuracy  of  thermionic  emission 
models.  As  similarly  detailed  experiments  are  not  available, 
we  do  this  analysis  with  our  calculations  for  RRKM-type, 
VW.  and  VW  +  IVR  models. 

We  first  consider  specific  rate  constants  kdel{E)  and  we 
ask  whether  kdc[{E)  and  fcdet,ce(T)  approached  each  other  when 
E  and  T  are  related  through  Eq.  (1,1).  kde{  x(T),  i.e.,  the  high 
pressure  limiting  rate  constant  of  the  falloff  curve  of  detach¬ 
ment  and  attachment1^  has  been  measured  in  Ref.  19  and 
exploited  to  derive  the  improved  value  of  the  EA 
=  1.20(±0.05)  eV.  WW,  kdel(E),  and  attachment  rate 
constants  and  cross  sections  in  Ref.  1 5  were  determined  in  an 
internally  consistent  manner  For  the  optimum  experimental 
temperature,  kdel  *(650  K)=  1 15  s”1  was  measured  in  Ref. 
19.  The  average  energy,  (£(650  K)),  at  this  temperature  is 
approximately  0.4  eV.  At  this  energy  obviously  one  has 
&det(£=0.4  eV)=0.  [On  the  other  hand,  &det(£=EA 
+  0.4  eV)*=  106  s“*.]  A  microcanonical  internal  energy  distri¬ 
bution  just  does  not  provide  energies  which  are  high  enough 
to  overcome  the  threshold  energy  EA=1.2eV.  The  finite 
heat  bath  concept  of  Eq.  (11)  in  this  respect  does  not  apply 
to  detachment  rates. 

We  next  consider  parent  reactants  SF6  *  directly  at  the 


moment  before  electron  emission  takes  place.  The  energy  E 
then  already  has  been  used  to  overcome  EA  and  the  remain¬ 
der  £-EA  is  partitioned  between  the  daughter  neutral  SF6 
and  the  electrons.  The  energy  distribution  /(£,e)  calculated 
in  Sec.  Ill  through  energy  conservation  is  complementary  to 
the  energy  distribution  of  the  daughter  neutral,  and,  on  aver¬ 
age,  one  has 

E-  EA  =  (Ed)  +  (e),  (4.1) 

where  ( Ed )  is  the  average  internal  energy  of  the  daughter 
neutral.  Identifying  ( Ed )  with  the  average  thermal  internal 
energy  of  a  canonical  distribution  of  SF6,  Eq.  (4.1)  allows 
one  to  define  an  effective  temperature  Td(E)  of  the  daughter 
neutral.  (It  is  common  practice  to  neglect  (e)  against  ( Ed ).) 
Then,  Td  follows  from  Eq.  (4.1)  with  (Ed)  obtained  from  the 
thermodynamic  relation 

5 

(Ed)  =  X  h  Vjiexpih  VjIkTj)  -  1  ,  (4.2) 

/=  I 

with  the  s=\5  frequencies  v}  of  SF6.  Often,  however,  Td  is 
not  derived  from  Eq.  (4.2),  hut  from  the  daughter  vibrational 
density  of  states  pd{E')  through 

kTd  =  [d\np/E')/dE']-'.  (4.3) 

This  simplification  introduces  inaccuracies,  the  latter  because 
pd(E’),  such  as  obtained  from  Beyer-S  wine  hart  counting,4  is 
an  irregular  step  function  over  the  energy  range  of  interest 
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FIG.  5.  Average  energy  (e)  of  electrons  detached  from  SF6~  at  lolal  energy  E,  Full  curves  =  detailed  calculations;  (a)  RRKM-type  model;  (b)  electron  capture 
VW-model;  (c)  electron  capture  with  IVR  =  VW  +  IVR  model  Dashed  lines= thermionic  emission  results;  (d)  pulling  (e)=0  in  Eq.  (4.1 );  (e)  including  finite  (e) 
from  Eqs.  (4.1),  (4.2),  and  (4.4);  see  text. 


and  needs  appropriate  and  grain-dependent  smoothing.  We 
have  avoided  this  practical  problem  by  using  the  Whitten- 
Rabinovitch  approximation  of  Eqs.  13.7),  (3.8a),  (3.8b), 
(3.8c),  and  (3.8d),  see  above. 

Td(E)  is  the  central  quantity  of  the  further  analysis  dis¬ 
cussed  here.  One  may  ask  whether  finite  heat  bath  arguments 
apply  and  allow  one  to  obtain  ( e(E ))  from  the  relation 

<«(£)>  ~*rrf(£).  (4.4) 


(which,  however,  may  not  be  realistic  for  larger  energies,  see 
above). 

Td(E)  has  also  been  used  to  express  k^(E)  in  the  sim¬ 
plified  form.6 

kJtE)  -  K£)Pj(£-EA)/P/,(£),  (4.5) 

with  densities  of  states  p,y(£')  and  pp{E' )  of  the  daughter 
neutral  and  the  parent  anion,  respectively.  v{E)  here  is  given 
by 


In  order  to  answer  this  question,  Fig.  5  compares  (e(£))  from 
Eqs.  (4  1),  (4.2),  and  (4.4)  with  our  detailed  modeling  re¬ 
sults.  One  realizes  that  Eq.  (4.4)  with  kTd{E)  from  Eqs.  (4.1) 
and  (4.2)  [i.e.,  without  the  common  assumption  (*-;)= 0  and 
without  using  Eq.  (4.3)]  gives  nearly  perfect  agreement  with 
the  detailed  results  for  the  RRKM-type  model  (curves  a  and 
e  in  Fig.  5).  On  the  other  hand,  neglecting  (*:)  in  Eq.  (4.1) 
leads  to  curve  d  in  Fig.  5  which  is  between  the  results  for  the 
RRKM-type  and  VW  models.  The  differences  to  the  VW 
+  1VR  results  (curve  c  in  Fig.  5)  are  even  larger.  Although  (e) 
from  the  finite  heat  bath  calculations  with  the  purely  thermo¬ 
dynamic  value  of  (ESf  (Td))  gives  semiquantitative  agree¬ 
ment  with  the  detailed  modeling  results  and  certainly  can  be 
used  for  simple  estimates,  dynamical  contributions  to  elec¬ 
tron  emission  such  as  specified  by  the  RRKM-type,  VW, 
and  VW  +  IVR  calculations,  lead  to  deviations  from  Eq. 
(4  4).  These  are  most  pronounced  for  the  VW  +  IVR  model 


v{E)  =  2/x[kT,,(  E)]2ctc(E)/  irfi  ■’  (4.6) 

and 

(tc(E)  =  I  £o-al(e)exp[“  elkTd(E)]cie 
Jo 

\  e  exp[-  elkTd(E)]def  (4.7) 

-M) 

with  the  attachment  cross  section  <ral(£).  We  have  explicitly 
calculated  v(E)  using  Td(E)  and  exploiting  the  relation 

E(TjE)  =  P(E,E)[kTd(E)]2l7rf},  (4.8) 

see  Secs.  II  and  III  Using  the  (slightly  corrected,  see  above) 
Whitten-Rabinovitch  approximation  for  pJLE')  and  p  (£')  in 
Eq  (4.5)  then  leads  to  the  results  which  in  Fig.  6  are  com¬ 
pared  with  the  original  kdcl(E).  There  are  differences  of  about 
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FIG.  6.  Specific  rate  constants  /^(E)  for  electron  detachment  from  SF6"  [full  curves:  Detailed  calculations  identified  in  Fie.  2;  dashed  curves:  Simplified 
thermionic  emission  approximations  lo  the  upper  three  detailed  cases:  The  highest  dashed  curve  is  from  Eqs.  (4.5)— (4  8),  the  next  uses  (slightly  corrected) 
Whitten-Rabinovitch  densities  of  states  for  pd,  and  the  lowest  curve  uses  Td  from  Eqs  (4.1),  (4  2),  and  (4  4),  see  text]. 


a  factor  of  1.5,  regardless  which  of  the  RRKM-type,  VW,  or 
VW+IVR  models  for  P(e)  is  preferred.  As  Eq.  (3  9)  gives 
better  accuracy  and  is  simpler  to  use  than  the  detour  via  Eqs. 
(4.5)-(4.7)  and  employing  Td{E),  we  do  not  see  advantages 
for  using  Eq.  (4.5)  for  the  present  system.  Instead  of  Eq. 
(4.8),  often  experimental  nondissociative  attachment  cross 
sections  have  been  used  in  Eq.  (4.7).  In  the  presence  of  VEX, 
then  detailed  balancing  may  have  been  violated,  see  above. 

Apart  from  Td(E)  there  is  still  another  effective  tempera¬ 
ture  in  use,  an  “isokinetic  electron  emission  temperature” 
Te(E)  This  is  defined  through  the  equation 

**,(£)  =  K£)exp[-  EA lkTe(E)].  (4.9) 

Having  determined  kdcl(E)  previously  and  taking  v(E)  from 
Eqs.  (4.6)  and  (4.7),  one  derives  Te(E)  from  Eq.  (4.9).  The 
results  are  shown  in  Fig.  7  together  with  7^(£)  based  on  Eq. 
(4.1).  The  effective  temperatures  Td{E)  and  Te{E)  defined 
above  in  Ref.  21  have  been  calculated  for  the  present  system 
as  well.  The  results  for  Td(E)  on  the  basis  of  Eq.  (4.3)  mark¬ 
edly  differ  from  the  results  shown  in  Fig.  7  (about  a  factor  of 
2  at  £-EA=0.05  eV,  a  factor  of  0.8  at  £-EA=0.5  eV).  In 
contrast  to  this,  the  results  for  Te{E)  (for  EA=1.2eV)  are 
accidentally  in  reasonable  agreement. 

If  pd(E')  and  f>p(E’)  were  similar  functions,  E  and  EA 
were  not  too  distant,  and  the  heat  capacities  Cd  defined  by 

Cd  =  d(Ed(Td))/dTd  (4.10) 

were  constant,  then  an  approximate  relation  between  Te{E) 
and  Td(E)  could  be  derived.  According  to  Refs.  6  and  9,  one 
would  obtain 


T'e.appCE)  *  T/E)  -  EM2CJE)  -  EA2l\2C2d(E)T/E). 

(4.11) 

The  evaluation  of  Eq.  (4.11)  for  the  present  system,  however, 
gives  incorrect  results.  For  example,  Te  app(£)/7f(£)  =  0.41, 
0.049,  -0.28,  -0.79,  and  -2.15  would  be  obtained  for  E 
-EA=  1.11, 0.8 1 ,  0.63,  0.44,  and  0.20  eV,  respectively.  This 
failure  of  Eq.  (4.11)  is  easily  understood,  none  of  the  condi¬ 
tions  used  in  its  derivation  are  valid  in  the  present  case.  For 
example,  the  difference  of  the  frequencies  of  SF6  and  SF6“ 
introduces  an  error  of  a  factor  of  580  into  the  ratio  of  pd  and 
pp  and  the  specific  heat  Cp  is  not  at  all  independent  of  Td(E). 
The  effective  temperature  Te(E),  therefore,  does  not  appear 
to  be  a  useful  quantity  for  the  present  example. 

The  transition  from  sparse  low-  to  quasicontinuous  high- 
energy  behavior  should  be  expected  to  take  place  when  the 
energy  of  the  daughter  neutral  exceeds  its  zero  point  energy. 
For  SF6,  the  latter  amounts  to  £,=0.583  eV.  In  reality,  how¬ 
ever,  the  smoothing  effect  of  electron  capture  theory  in  com¬ 
parison  to  RRKM-type  calculations  of  fcdet(£)  reduces  the 
maximum  deviations  from  the  smoothed  curve  to  ±2%  al¬ 
ready  at  £-EA~£,/2,  see  Fig.  2.  In  contrast  to  this,  the 
maximum  fluctuations  of  f)d(E=Ez)  still  amount  to  more  than 
±50%.  This  is  reflected  in  large  fluctuations  of  A £(£,e,)/Ae 
such  as  illustrated  (for  smaller  energies)  in  Fig.  1(d).  In  any 
case,  transitions  between  low-energy  and  high-energy  behav¬ 
ior  take  place  when  £-EA  is  roughly  of  the  order  of  £,. 

We  finally  comment  on  the  difference  between  the  val¬ 
ues  for  EA(SF6),  obtained  from  the  third  law  analysis  of 
measured  thermal  detachment  and  attachment  rate 
constants1^  (EA=1.2eV)  and  from  a  fit  to  measured  SF6" 
lifetime  distributions  on  the  basis  of  thermionic  emission 
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FIG.  7.  Comparison  of  the  daughter  lemperature  TJ^E),  based  on  Eq.  (4.1)  and  an  effective  or  isokinetic  temperature  Te(E ),  from  Eq.  (4.9).  in  thermionic 
emission  models  of  electron  detachmeni  from  SF6”  (delailed  calculations  of  this  work,  see  text). 


expressions  (EA  =  1.4  eV).  The  former  approach  employed 
the  experimental  attachment  cross  sections  from  Refs.  24  and 
25  such  that  a  comparison  with  £de,(E)  from  the  VW  +  IVR 
+  VEX  model  should  be  made,  see  the  lowest  curve  in  Fig.  2. 
Employing  EA=1.2eV,  in  Ref.  21  a  factor  of  five  times 
higher  values  for  £de,(E)  were  derived  than  in  the  present 
work.  A  factor  of  2  of  this  difference  is  due  to  the  omission 
of  one  of  the  two  electronic  degeneracy  factors  [see  Eq. 
(2.1 1)  versus  Eq.  (4)  from  Ref.  21].  Another  factor  of  about 
1,5  is  the  consequence  of  using  the  simplified  Eq.  (4.5)  in¬ 
stead  of  the  original  Eq.  (2.1)  for  fcde,(£).  This  is  illustrated 
by  the  differences  of  detailed  and  simplified  results  in  Fig.  6. 
Furthermore,  neither  the  experimental  cross  sections  nor  the 
VEX  factor  of  Eq  (3.4)  is  expected  to  be  valid  over  large 
energy  ranges,  apart  from  the  fact  that  experimentally  ob¬ 
served  attachment  cross  sections  and  kdc]{E)  most  probably 
do  not  correspond  to  reverse  rate  processes.  Experimental 
attachment  cross  sections  and  calculations  with  the  VW 
+  IVR  +  VEX  model  then  both  should  not  be  used  for  calcu¬ 
lations  of  *dc«<£>. 

V.  CONCLUSIONS 

The  present  application  of  thermionic  electron  emission 
concepts  to  the  specific  rate  constants  kdel(E)  and  energy  dis¬ 
tributions  /(£,/-:)  of  the  detachment  of  electrons  from  SF6“ 
have  provided  insight  into  the  accuracy  of  the  currently  used 


models  It  is  emphasized  that  attachment  and  detachment  are 
only  linked  by  detailed  balancing  when  they  correspond  to 
truly  reverse  processes.  The  comparison  with  calculations 
from  statistical  rate  calculations,  combined  with  electron 
capture  theory  and  accounting  for  kinetic  complications,  re¬ 
veals  considerable  inadequacies  of  the  simplifications  used 
within  current  thermionic  emission  models,  at  least  in  their 
application  to  the  considered  SFft”  system. 

The  present  work  demonstrates  that  the  detailed  state- 
resolved  calculation  of  £del(£)  from  Ref.  15  can  considerably 
be  simplified  when  (slightly  corrected)  Whitten-Rabinovitch 
expressions  for  the  densities  of  states  are  used,  see  Eq.  (3.9). 
However,  the  gain  in  simplicity  is  at  the  expense  of  losing 
the  low-energy  fine  structure  in  £de,(£).  It  is  illustrated  that 
continuous  energy  distributions  of  the  form  of  Eq.  (3.10), 
i.e., 

/(E,e)  x  (£/<£»"  exp(-  £/<£» 

are  obtained  which  have  a  “nonthermal"  shape,  i.e.,  n 
=  0.15  is  found  which  is  markedly  smaller  than  “thermal 
expectations”  of  n  =  0.5  or  1.  The  average  energy  (*;)  as  a 
function  of  the  total  energy  E  also  shows  dynamical  (and 
kinetic)  influences  of  the  process. 

Apart  from  the  dynamical  refinements,  the  relation  be¬ 
tween  (e)  and  E  in  accord  with  the  thermionic  emission  con¬ 
cept  can  be  roughly  estimated  by  Eqs.  (4.1)  and  (4.4),  i.e.,  by 
the  relation 
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£-EA~<£SF6(7;)>  +  <e>,  (5.1) 

where  ( £SF  (Td))  denotes  the  internal  vibrational  energy  of 
SF6  at  a  canonical  temperature  Td  and  (t')~kTd.  The  depen¬ 
dence  of  Td  on  E  can  easily  be  calculated  through  Eq.  (5.1) 
which  corrects  the  sometimes  used  practice  of  omitting  the 
last  term  kTd.  The  representation  of  fcde1(£)  in  terms  of  Td(E) 
appears  less  useful,  as  the  direct  calculation  of  *de.(£) 
through  Eq.  (3.9)  is  easily  done  and  avoids  unnecessary  sim¬ 
plifications  which  introduce  errors.  Equation  (3.9)  also  as¬ 
sures  the  intrinsically  consistent  detailed  balancing  link  to 
attachment  properties.  At  least  for  the  present  system,  the 
introduction  of  an  “isokinetic”  electron  emission  temperature 
Te(E)  appears  less  useful  because  the  simplified  relation 
(4.11)  between  Te(E)  and  Td(E)  in  the  present  case  was 
shown  to  lead  to  erroneous  results.  The  effective  temperature 
Te(E),  therefore,  does  not  appear  to  provide  much  insight 
into  the  details  of  the  phenomenon.  The  temperature  Td(E) 
from  Eq.  (5  1)  arising  from  the  thermionic  emission  concept 
instead  appears  more  useful,  although  it  misses  the  dynami¬ 
cal  details  of  the  processes.  Nevertheless,  it  provides  a  direct, 
yet  approximate,  estimate  of  the  average  electron  energy  (e). 
One  should  note,  however,  that  the  shape  of  the  relative  elec¬ 
tron  energy  distribution  f(E,e/(e))  markedly  differs  from 
expectations  for  thermal  distributions. 
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APPENDIX:  MOLECULAR  PARAMETERS 

The  calculations  of  this  work  employed  the  vibrational 
quanta  (in  cm”1)  for  SF6  and  SF6~  from  Ref.  22.  They  are 
given  by  346(3),  519(3),  611(3),  655(2),  779(1),  and  965(3) 
for  SF6,  as  well  as  237(3),  336(3),  306(3),  447(2),  626(1), 
and  722(3)  for  SF6”,  with  degeneracies  shown  in  parenthe¬ 
ses. 
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